A non-local OPE for hard QCD processes near the elastic limit 
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A leading twist expansion in terms of bilocal operators is proposed for the structure functions 
of deeply inelastic scattering near the elastic limit x — * 1, which is also applicable to a range of 
other hard quasi-elastic processes. Operators of increasing dimensions contribute to logarithmically 
enhanced terms, which are suppressed by corresponding powers of 1 — x. For the longitudinal 
structure function in moment (N) space all the logarithmic contributions of order ln fe N/N are shown 
to be resummable in terms of the anomalous dimension of the leading operator in the expansion. 

I. INTRODUCTION 

The treatment of hard processes in QCD perturbation theory involves the factorization of the short from the long 
distance dynamical regimes, characterized by mass scales Q and Aqcd respectively. Such processes become quasi- 
clastic in those regions of phase space where a final state invariant mass M becomes much smaller than Q. For the 
perturbative treatment of these three-scale processes the condition M 2 ^> Aq CD is also necessary. Typical quasi-elastic 
processes are deeply inelastic scattering (DIS) at large Bjorken x, Drell-Yan production near the partonic threshold 
and thrust in e + -e~ annihilation near the two-jet limit. We argue that existing methods for studying quasielastic 
processes can be extended in a systematic way, via a non-local operator expansion to be introduced below. The DIS 
longitudinal structure function will be considered for dcfinitcncss. 

The longitudinal structure function (Ft,), like any DIS observable at large Q, may be treated by the Wilson 
OPE, organized through the light cone expansion or, equivalently, factorization formulas. Given regulator (e) and 
factorization scheme (MS), the leading twist factorization at scale fj, 2 can be written as 



F L (x, Q 2 ,e) = C L {x, Q 2 //i 2 ,a s (/i 2 )) ® f(x, a s (^ 2 ), e) + 0(A 2 /Q 2 ) , 



where we define the convolution in longitudinal momentum as 



h(x) <g> f(x) = / dxi / dx2 S(x — xix^) h(x{) /(X2) , 



(1) 



(2) 



and the perturbative expansion of any function R as R — 'Y^ =f) {cy. s /'n:) n R^ n \ The coefficient functions can be 
computed in perturbation theory and the results to the second order in the coupling are |l| |3[ 



4 0) (x,Q 2 / M 2 ) = 0, 



C l 1) (x,Q 2 /ii 2 )=C f x, 



C L 2 \x,Q 2 / i ?) = 



In- 



(3) 



with P qq the quark splitting function and (3\ the one-loop coefficient of the beta function. The coefficient function at 
Q 2 = [i 2 is 



C[ 2 \x,l) = ^-C 2 F \n\l-x) + 



1 



2 l„2. 



4 - W 



C 2 F ln(l - x) 



(C(2) - l)C F C A - X -fi x C F 



ln(l — x) + rcg 



(4) 



Only terms non-analytic at x — 1 are displayed in the above equation. Our objective is to show that these terms 
can be captured by expectations of appropriately defined operators, and that they can be resummed to all orders in 
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perturbation theory. It is instructive at this point to compare the above coefficient function with the corresponding 
one for F 2 to 0(a s ), 

CP (x, Q 2 /^ 2 ) = \P$ In ^ + 4C F i^Y^) + ~ 3 (1^)7 + ^ ln (! ~ x ) + re S- ( 5 ) 

We note the absence of plus distributions in the case of Fl . The longitudinal structure function has only logarithmic 
divergences to all orders in perturbation theory. Upon taking the Mellin transform 

F(N,Q 2 )= [ dxx*- 1 F(x,Q 2 ), (6) 
Jo 

the leading logarithms of Fl are of order ln k N/N in moment space. Near the elastic limit (N — > 00) they are power 
suppressed relative to the order ln fc N terms coming from the plus distributions of F 2 . For this reason, only recently 
has the Sudakov factorization formalism been extended to capture purely logarithmic threshold corrections 0] . 



II. THE NON-LOCAL OPE IN THE ELASTIC LIMIT 



The Sudakov factorization formalism of perturbative QCD applies to hard processes in the quasielastic regime Q 2 ^> 
M 2 3> Aq CD and resums terms of order In N, where N is the moment with respect to the variable (Q 2 — M 2 )/Q 2 . 
We shall refer to this formalism as the leading-jet approximation, and we briefly summarize it here before extending 
it to the case of Fl- (For a recent review and references see Ref. B). 

Every DIS observable is obtained as a particular projection of the hadronic tensor 

(P. ?) = ^ / d 4 y e- 1 ^ (p\jl (0) j v (y) \p) . (7) 
Structure functions are obtained as Fl = P^Wuv, F2 — P 2 V Wu,v with the projectors 

p r = ^W, ir=-*r+f*r- («) 

The leading-jet expansion of F 2 leads to the following factorization formula 

F 2 (x,Q 2 ) = \H 2 {Q 2 )\ 2 j dx' dw J((x' ~x-w)Q 2 ) V(w) cj)(x') (l + 0{l~x)Y (9) 

All three factors in the longitudinal momentum convolution integral are defined as expectation values of non-local 
operators. Explicitly, the soft function is 

V(w) = J dy-e-*"*'** (0|$t (0, -00) $ c (0, y~v) -oo)|0) , (10) 

the initial state jet function is 

<P(x') = J dy~ e-"'^- s (7#(0) jb $ c (0, y~) j(y)\p) ® V' 1 , (11) 

and the final state jet function is 

J ((1 - z)Q 2 ) = IdHe-^^ v <0|$t(0, -00) V(0) $(y) $ v (y, -00) |0) <8> V^ 1 . (12) 



We have defined the = light-cone vector in the direction of motion of the incoming parton and w M as the opposite 
(— ) direction. We have also introduced the Wilson line operator 
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<J>. u (x + tv, — oo) = Pexp 



ig s / ds A fl (x + st>) 



DO 



(13) 



The perturbative expansion of the Wilson line operator generates the Feynman rules for the eikonal interaction 
between a jet of fast moving partons in the v^ 1 direction and the soft gluons. The presence of the V~ x factor in the 
jet definitions removes any purely soft contributions from these functions. In a perturbative calculation of the jets, 
this corresponds to the prescription of applying soft subtractions in the contributing diagrams order by order. 

In this operator language we can see that the factorization formula Eq. (^) comes from a rearrangement of the 
fields that define the currents to which the DIS probe couples (in our case the electromagnetic current j M = ipj^i/j). 
Indeed, by "wiring" one of the fermion fields to the initial state and the other to the final state we obtain 

PfWp, = PrFT( 4 )(p|^(0) 7 ^(0)V;(2/)7^(j;)|p) « FT^(p\^(0)^(y-)\p) ® FT^F^ W(0)^2/)|0) 7 , , (14) 

where we have introduced FTii? and FT^ 4 ^ as shorthands for the one- and four-dimensional Fourier transforms. In 
this, leading jet, approximation, the F r , r = 2, L, DIS structure functions are written as convolutions of two jets, which 
are identified as the initial and final state jets defined in Eqs (|ll],[l^). Note that this procedure would give a vanishing 
result for Fl, because the projection of the final state matrix element with Pl, Eq. (||), is zero by the Dirac equation 
of motion for the on-shell incoming quark. Fl vanishes in the leading-jet approximation, and the coefficient functions 
of Fl are less singular in the small parameter 1 — x relative to the ones of F%, as we noted in the introduction. The 
operator formalism indicates how it can be expanded to capture such terms that are suppressed by powers of 1 — x. 
To this end we introduce higher dimensional operators with the quantum numbers of the electromagnetic current and 
rearrange their fields between initial and final states. Let Of (y) = "0(y)Sf (y) be an operator of dimension higher 
than 3. Then, as in Eq. (jl4|), we obtain a new contribution 

prw^ = pr \h l (q 2 )\ 2 FT^H P \4,m»Ao)ziAyMy)\p) 

« \h l (q 2 )\ 2 FT«<2#(oh%-)b> ^FT^ipPrm^mlM ) ■ ( l5 ) 

with Hl{Q 2 ) a new coefficient function. The second factor in the convolution product above is a new jet function. At 
each mass dimension for the operator Of we generate a class of jet functions whose contributions are suppressed by 
one factor of 1 — x relative to those of lower dimension level. Therefore, an expansion in powers of 1 — x reduces to 
dimensional counting, subject to the equations of motion of the incoming and outgoing partons. This is the non-local 
(bilocal) operator product expansion near the elastic limit ||. In the language of effective field theories, we have 
"integrated out" the dynamics between the scales (1 — x)Q 2 and Q 2 , which generates in a new operators Oi in an 
effective Hamiltonian for the still-perturbative dynamics at the intermediate scale. 

The bilocal OPE in the elastic limit is a method of capturing corrections of order (1 — x) n in the final state 
interactions. This expansion is still leading twist in Q 2 . This may be understood by observing that, if the dimension 
of the effective operator Of is 3 + n, then a field of dimension 3/2 + n is emitted into the final state and the phase 
space integral scales as (1 — x) n Q 2n . It is therefore the phase space integral that absorbs the factor 1/Q 2n associated 
with the coefficient functions for the higher dimensional operators, leaving the factor (1 — x) n , which generates an 
expansion in 1 — x. 

Let us apply the above procedure to Fl- We seek to define a new jet function that captures the leading power 
contributions in 1 — x. The dimension 4 operators with the quantum numbers of the current are 

0% a = v*$p ± 1>, 0& = v"#xV, 0£ C = ^D W <>, 0% d = $D^. (16) 

It is easily seen that On, and 02 C have vanishing projection when contracted with P^ . They do not contribute to Fl 
but can contribute to F2 at the subleading level. Oid has a longitudinal projection that is proportional to ip{v- D)ip- 
After integration by parts this term vanishes by the equation of motion of ip at x — 1. We therefore conclude that 
the only operator of dimension 4 that will contribute to Fl at the leading level is 02 a , with longitudinal projection 



0$ a = $p ± i>. (17) 
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This results in the jet function for Fl 



J'((l - z)Q 2 ) = (J- ^\ FT$ zp (0|$t (0, -oo) Hv) f> ± -oo)|0) ® V" 1 , (18) 

where we have restored the Wilson lines for manifest gauge invariance and have divided out any purely soft contribu- 
tions. The resulting factorization formula for Fl is similar to the one for F 2 . 



F L (x,Q 2 ) = \H 2 (Q 2 )\ 2 J dx' dw J' ((x' - x - w)Q 2 ) V(w) <p(x') (l + 0(l-a;)) 



(19) 



In the Fl resummation formula the anomalous dimension of J' will appear. 



III. BEYOND THE EIKONAL APPROXIMATION - RESUMMATION 

The factorization formula for Fl at large Bjorken x was derived via the non-local OPE near the elastic limit in 
the previous section. It is however instructive also to study the problem from the point of view of perturbation 
theory diagrams. This approach yields the Feynman rules for computing the anomalous dimensions of the higher jet 
functions. 

The starting point of the diagrammatic approach is the observation that both F 2 and Fl are different projections of 
the same hadronic tensor . Any non-analytic terms in their expansion in massless perturbation theory arise from 
well known configurations in loop momentum space, the Landau pinch singular surfaces of W^ v . We can therefore use 
the infrared power counting and the jet-soft factorization techniques on which general QCD factorization theorems 
are based (for a review see Ref. 0). Given the light-cone direction vectors and v M , normalized as v • v = 1, the 
momenta of the incoming quark and the outgoing jet p and p respectively are 

'-v-'fa + fa- (20) 




pf 

FIG. 1. The factorization of the u-collinear gluon (I) from the hard subgraph generates the effective vertex Cg>. Parallel lines 
around the photon denote contraction with p M . 

Consider the configuration shown in Fig. 1, in which the depicted gluon belongs to the v-jet. The fermion line 
carrying momentum p — I is far off-shell and belongs to the hard scattering subdiagram Hl- If, however, we were to 
factorize the top collinear gluon by eikonalizing the fermion propagator, then this, in general leading, IR contribution 
would vanish in Fl because of the contraction with p M from at the current vertex. Note also in the case where 
the gluon is soft, the leading eikonal contribution is similarly projected out. Therefore, the expansion will start with 
terms that are subleading in the usual eikonal approximation. Expanding the diagram in the gluon's virtuality we 
obtain 

9s^ TTT2 m~ f)l X - 9s ( 7 i - I, ■ 7X^) - OHD , (21) 

where we have dropped terms vanishing by the incoming quark's equation of motion. An effective vertex to 0(g s ) 
has emerged. This vertex, unlike the usual eikonal ones, probes the transverse momentum of the struck quark. It 
satisfies the following algebaic properties, that are useful in diagrammatic calculations: 
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O x (l)h = 0, O x (l)v x = 0, O x (l)v x = -g s O x {l)O x {l') = -2g 2 s , O x (l)...0 X (l') = sM-TiA • (22) 

v ■ i 

Moreover, it is noted that the vertex O x is the 0{g s ) expansion in momentum space of the operator $„(0, — oo)p±ip(0), 
i.e. half of the bilocal operator that defines J', Eq. (|Tg|). It can be shown by induction that the above construction 
generalizes to any number of collinear gluons (?]] . The factorization of the final state collinear partons from the initial 
state jet defines a final state jet function exactly as in Eq. (^8|). 

Resummation may be thought of a consequence of factorization. |^] Once an observable is written as a product of 
effective operators, then the UV renormalization of these operators generates the evolution of the coefficient functions 
of the observable. Convolutions become simple products in moment space where the Mellin transform of Fj,, Eq. (|^), 
can be written as 

/ ( I - \2 \ 2 



(p ■ v) 2 (p ■ v) 2 2 , 
Hl I - — , - — , a s (fj, , 



i ~ v ( Q 2 (P-vf 



X—J 

N 



Nfi 2 



ViV 2 /x 2 



(23) 



<iz z 



N-lj, 



(i-z)g 2 ( P -.«) 5 



,a 5 (/J 2 



(24) 



Here, the factor of 1/N has been exhibited explicitly, 
1_ ~, (p-.f 2 N r 1 

N J {Nfx 2 ' M 2 ><W )j Jq v ^ ^ 

The treatment of the J' evolution proceeds in the usual way ||. The jet function satisfies two equations. The first is 
the Lorentz covariance equation, which describes how the jet varies in response to changes in the light-cone direction 
with respect to which it has been defined, and the second is the jet's own renormalization group equation, 
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d ln(p • v) 2 



In J 1 



Q 2 (P-v) 2 ( 2 A K ( Q 2 , 2 



N/i 2 
d 

din /j, 2 



\N^ 2 



(V('^^,«V))- 



~ ( Q 2 (p-v) 2 2 \ l f( 2^ 

In J I ~rz ~ i — ,a a (fi ) =--7J'(cW ))• 



\N/j, 2 



(25) 



(26) 



The right-hand side of the covariance equation is renormalization group invariant and the K factor renormalizes as 

(27) 



where is the well known cusp anomalous dimension in MS scheme, 



IK 



7T V 7T / 



0(al 



Eqs (|25| , p6| ) can be solved simultaneously subject to the initial condition 

a s (Q 2 "> 



J'(l,l,a s (Q 2 )) = C F - 



0(ai(Q 2 )). 



(28) 



(29) 



from the Mellin transform, Eq. (24), of J' to lowest order 0{a s ). To this order, there is only one contributing diagram, 
depicted in Fig. 2. 




FIG. 2. The lowest order contribution to the jet function J'. 
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Substituting the solution of the jet equations (|25|) and (|26|) into the factorization formula (f23|), we obtain the 
resummed expression 



F L (N,Q 2 ,e) = C F ^P-(V4>)(V 2 =Q 2 /N,e) 



x exp 



ttN 

>h ^ln N 1K (a s (m 2 ) ) + TP (a, (fi 2 ) ) + 2(3(a s ) ^- In J' ( 1 , 1 , a s ( M 2 ) ) 



>/N A 4 







ln°jV 

AT 



(30) 



Apart from the jet J', A~ dependence is also found from in the factor V</> in Eq. (p3|) . This dependence, however, 
cancels in the coefficient function Cl once we divide Fl by the Mellin transform of the parton distribution function 
/, Eq. (Q). In Eq. @, the contributions coming from the UV renormalization of the jet function J' are contained in 
7 j' and are kept distinct from the renormalization of the coupling, contained in the /3-function term of the exponent. 
This distinction is important for the bookkeeping of the renormalization terms in the perturbative calculation of the 
anomalous dimensions. 

In summary, we have shown that the extension of the bilocal OPE to the case of Fl modifies the Sudakov resum- 
mation formula essentially only through the introduction of the anomalous dimension of the new jet function J'. The 
cusp anomalous dimension contains all universal soft radiation effects. 



IV. THE JET ANOMALOUS DIMENSION TO O(ALPHAs) 

Inspecting the resummation formula, Eq. (^o|), we see that the only new ingredient is the anomalous dimension jjr. 
In this section, we outline its calculation to 0(a s ). Because Fl itself and J' start at 0(a s ), F^ = J'^ = = Cfx, 
we need to compute J' to 0(a 2 s ). Using the rules for the effective operator O x of the previous section we identify the 
second order diagrams shown in Fig. 3. 




FIG. 3. The 0(o? s ) diagrams for the jet function J'. 

First we compute both the UV and IR singular parts of the diagrams. Then we perform an IR subtraction as 
dictated by Eq. JTs|). The infrared part to be subtracted out is J'^ <X> and we note that the soft function V is 



G 

i 



the same for Fl and F2. After the subtraction, only UV simple poles survive and from these we construct as usual 
the anomalous dimension 7,/' . The result is 



U'{Us) 



Q-s 
n 



-C F - 2C A - 4C(2) C F 



C a 



0(al 



(31) 



Comparison of the resummation with the fixed order results listed in the introduction is made by first writing the 
resummation prediction for the coefficient function to O{or s ) 



PV(2) _ C F 



C 



2N 



(1)1 2 
Ik ln 



N 



1 



N 



(32) 



with No — e 1B . Upon substituting for the anomalous dimensions and 77' from Eqs. (|2^, |T|) we obtain the Mellin 
transform of the fixed order result in Eq. (^). 



V. SUMMARY 

We have shown how the light-cone expansion for hard QCD processes near the elastic limit can be extended into 
a general bifocal OPE that captures power terms suppressed in 1 — x. This operator expansion is based on purely 
dimensional arguments, but it can also be derived from the analysis of the perturbative diagrams in the infrared 
limit. The result is a Sudakov resummation formula for the observable, in which the anomalous dimensions of the 
higher dimensional jet functions enter the exponent. We have analyzed the DIS longitudinal structure function at 
large Bjorken a; as a concrete example of this operator expansion near the elastic limit. The anomalous dimension of 
the Fl final state jet was computed to first order in a s , and the resummation formula for the coefficient function was 
shown to agree with the fixed order calculations. 

We emphasize that our analysis is leading twist in Q 2 and is consistent with the usual light-cone expansion, ft 
resums terms that are of order ln fc N/N, k > 0. This approach can be extended to the case of F2, where such terms 
are present, but are subleading relative to order ln fc N, k > terms. Indeed we expect that the operator expansion 
presented here is relevant for all hard QCD processes in the quasi-elastic limit. 
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